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(54) A UNEAR INTERPCMJimG METHOD FOR COLOR 
SIGNALS IN A MEMORY 

(71) W^ .DAINIPPON SCREEN SEIZO KABUSHHtt KAISHA, a 
Japan»e Corporatum, of l-l Tcnjln-kitaittaaiii, Tcranouchi-flgani 4-dioine, Horikawa- 
dcd, KaxnlgTO-ku, Kyotl-shi, Japan, do hereby declare the invention, for which we 
Dray that a patent may be granted to us, and the mediod by which it is to be per- 
foxnied, to be partxculariy described in and by the following statement:— 

This invention relates to a linear mteipolating method for signals in a memory 
, which is used for color correction of picture signals in a reproducbg machine such 
as a color scanner, a color facsimile producer, or the like, in which color separation 
. picture images are produced by photo-^ectric f^^n^ng 

In convemim^ color pbotognpliic plate making color correction b often made 
fay. photx^phic maddng. Hontvcr, tius method has xnany defects, for example: 
limttadons of color conectfon, ability, necessity for many skflled engineers, uzueliable 
results of the color separation, irregular qualhy of finish, complexity, and tbt like. 
. In order to overcome these defect^ a color correction masldng m^od by an 
' 1$ electronic color sq^aration machme sudi as a color scanner has been developed and is 
nowadays mor e popular, Most of the color scannen- now used mjrfoy an analog 
^T 'wnputer system for the .color cotiection calculations so as to tecrease the calculation 
speed. 

This meth o d, however, has also defects such as the difficulty of the introduction 
of many kinds of ealculations because of tise restdcdon of calculation ability, tncvkable 
tStOB of ten^ferature drift and noue, multiptlicity of operational amplifiers and so 
forth as electric e lemen ts, inoonvenimoe of qieradon i&e. to many adfustments of 
potentmmeters and switchcSi and high mazuifectuziog cc^ 

If the analog computer system is simply replaced with a digital computer system, 
25 which has advantages such as a wide correction variable rahge .and convenieoce of 
operation, the calculation speed for the color correction decreases very much, and the 
processing ability is reduced. Accortlzngly, ditt system is not practicable. 
• Rooaatiy, a dlrea scanner has been developed for plate makipg m prfnthig, which 
pwforms color sepantioQ, color correction, omvefsiwi of scale of die reproduced 
30 fanag^ and halftone processing at the same time so as to meet die requirement for 
. mgh quality printfaig and rapid operation. In dds case, however, there is the defect 
that supplementary n^askmg or hand retouching dtat the color separatum canhot be 
apimed, as opposed to conventional color scarming whidi includes color ^Twratimi, 
color correction, conventian of scale of the reproduction ima^e, and halftone processing. 
35 / . ^ general, an original color picture is scanned by a color a^nn^ to obtain diree 
(red, greai, and hlue) color separati(m sigoals. These three color sqiaraticn signals 
are sent to a color opcruion arcm't; thetdsy finally ohtahiing recording « gnM s for 
density of printmg inks, such u <7an, magenta, yellow, and black. 

In/order to provide the most accurate possible color rcproducdoa, a combination 
40 of the amounts of cyan, magenta, and yellow inks (the Uack inl^ and so forth, are 
omitted for the sake of brevity of explanadon) is necepri^ A^^^iyf^ correspondmg 
to a combination of red, green, and blue color separation ^i griaiy 

Consequendy, for die purpose of color correctioii by selecting the combuxation 
of cyan, magenta, and ycflow values corresponding to the combination of red, green, 
45 and blue ^mlues,* the color-corrected combinations of cyan, magenta, and yeDow values 45 
coneqkonding to each com bi n a tion of red, green, and blue values are stored in a 
memory in advance, and then the color-corrected combination of cyan, magenta, and 
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ycUow vahitt is read imt by addressing the memory by the combination of red, gieen, 
and blue values conespondmg Aereta 

If each red, green, and blue range is divided into, for example, two hundred tone 
steps, altogeAtt 200* =8,000,000 comMoadons of cyan, magenta, and yeflow vahies 
mist be stored m the memoiy, which requires tfutt the memory have a larec capadtv.* 
This meana highotst; and thus is not practicable. o r j 

ihcrcfor^ in order to reduce the storage capacity required for the memorVs each 
aJ!^^J^ "^l^V^ ^"^^ f« enmple. aiztccn ton/ stow, 

Ite toe ^ge capacftjf requuemimt for 4e memory is reduced to a manamHe 
icveL On tte other hand, the tone steps become too rough, and the lack of cSmlt 
consistency becomes conspicuous, so that printrnt quality suffers. Therefore, m M 
ofiCi It 18 necessary, to mterpolatc mtcrmediate values properly between each two tone 

SKpS* . 

^sJ^iT^ invento ^ to provide a finear interpolatiag method for such 
sstgnals hi a memory ftee from ^ above motioned defects, which enables to 
mttno^ to ratcukte quicWy by using a shnple formula interpoktioo values, witihout 
IXidESL^ °^ of the mterpolatcd values b«weeii one imapE 

Acrading to the present invention Aere is provided a linear interpolating 
method for colmr signals in a memory of « picture lepiodudng machme,.comprisS 
stormg ap^nate vahies of color picture ouqmt signals conespondhig to certain 
stepped vak^ of cdor input signals hi the memow addressed in a three-dimensional 
fa^ and mtcrpok^^^ values of cdor output signals at pomts wMch are between 
saia valUM^ dmdmg up the cubic interpolation umt* of the memory whidi is * 
consmmed by a single 5tq> of each of the cdor input sknals into a rfurality of 
tetiahedra whose vertices are either vertices of the cubic unit; ccntcn & its faces, 
or Its com pom^ calculathig the cdor output siignal at tach vertex of the^ teo^ 
hcd» which is a center of a face of the cubic uniii if any, by averaging the values 
of the color output signal at the four vertices whidi are coiners of^ face. ^ 
at tte c«te poirt of tte cuWc unit by avcragmg the 

at afl eight of the vertices of the cubic unit, detetmmmg whidi of these te^ahd^ 

^opdated, and dcnvmg the mterpdated vafae at the inte^^cmlSto u 
wci^ Sim of the vahies at the four vertices of the deteimSSi tctrai^tf^ 
jfTJlAlS^iy^ ^1 oorrespondhig to the ratio dSrvda^ 
^ verdccs are the hitetpolatfon pomt and the other three 

vertices of the determmed t^rahednm to the vdume of the dl^ennmed tetrahedron. 

In one embodhnent in accordance widi Ae 'Svoaira the ciiic iioi 'ii'SiitiSd 
mto twenty-four tetiahedra, of each of which one vertex is the center of the cubic unfL 
two vcrfic« are vwtiocs of tbc cubic unit are ccamected by an edge of Ae cubic 
um^ and the fourth vertex is tte center of a square face of the diWc um't one edge 
of which face is die said edge dtiie cubic unk. «*w«*ec 

In a ^ccnd embodiment the cubic unit is dMded into six tettahedra by three 
planes which have a line m commoi^ said Unc being a long diagonal of the cubic unit, 
•and each plane containkg two eclgea, and four vertices of the cubic un^ 

In a ^tffl fimhtf embodhnent the cubic unit Is divided into five tetrahedta by 
four planes, each of whidi contams exactly three vertices of die cubic unit said 
pianes intersectu^ one another ahmg Imes which are diagonals of the fiices of die 
CQOzc unit. 

The present mvortfoa will nov be desoibed in detiil with i«^rt fn r*.. " 
accompanying drawings, in which: ««««« m aeiau wita respect to the 

^f-.i''^ \ of « convemioMl intopolating method over a two- 

dimensnmal mteipolatioa unit square; gvcr a two- 

o^JTJJS^I «re scheniatic Wewi of a squaie ioteipalation unit renbn and a 
SLi^aSSaSnK.'' .wosUmensiond SS" thre. 

tiond''Si5Lt^£2:2i,25^ val™, of convn. 

^^^6Uj^schen«tfcW 

7 i« a schematic view of a aiKc ■Dtetpohaian unit disected into twenty-four 
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tetrahoini aocoidirig to one of the miaiions of the method of the present invcation; 

Fjg. Jis aschenwtk vicwof one d.ihetetnibcdraof thediss^ of Fig. 7; 
. . ; 9 is a scheoiatk view of .another variatioa. of the method of the present 
iiweottoQy Wheick the unit cube fs dissected Dito su 
5 : vYFig."10 is* a schematic view of the tetiahcdia obtained by the dissection of Fig. 9; 

ana *'i • * ' ' 

' Figs. .11 and 12 iOustrHte anodicr method of dissecting the unit cube into -five 
tctjvhpdra, which gives anodier variation of the method of the present invention. 
/ In Qtder that the nxthod may be t>eiiet undexstoodi some cxplaiiatira of prior 
10 art.metfaods of tntmoladon will now be g i v en. 
. fi«fcfring to r, jhajt is dwOT 
of two variables, wtoe the mterval m be interpolated over is taken as unity. 
' Ai»r^ ^ y^¥L^ 3r*+yf) at a^point. P in aft interpolation region 

ABCD will be found by % mefhem a nca l interpolatipg fnediod, io which and yi are 
the mtegral pans of X aoul y and if and y, arc the decanal parts. 

For the imeipolation it is necessary tto the function at ttc vertices A» B, C 
m^f^-fr?^^^^ ^u'^tTV^Y'^i^J^'^i^jrt), U(ii+l.y,+l), ami 
V J}f mterpolated vajuc U(xor) wiU be a function of x^, y^ Uridyl), 

. U Xi+l^O, U(Xf+l^.+l), and U(«uyi+1). Further, for consii^cy, the ffi 
polatttl value should be consistent widi the known values of the original function at 
the QDtnexs of the unit region. 

>^ inteipolating method sadsifying such a condition will be described. It a 
<^llctt Imcar. mceipolation because on the edges of the unit region it reduces to a 
' simple Imear intexpolatian function. 

? ^ ^ value. XJ(j^) at the p&int F in the mtcrpolation unit uuaie 25 
' -I^i . ™* perpendiculars fcom the pomt P to each side AB, BC CD, 
A ^^Ll^ Designate the feet of these perpendiculars by Qi, <i, and 
ic^eaively» as shown in Fig. 2, and add up the results' obtained by multiplying 
ea^ known- vahic at the vertices A, B, Q and D Iqf the area of each rectangle opM 
tn the vefteac,the|«by obtaining the foHovring equation (^^ • 30 

V(w)=U(si+3^yiH-yf)=^UCibyk) (1-x^) (i-^^^^ 

HJ(Ki+ljrO XI (1-yO . (I) 
+U(ii^,+1) (1-xi) J, 
-HU(ii+lo^+l) xg y, 

The tnterpolattng method according to the formnia (I) satisfies the above 
boundary conditions at the comers of the unit square and reduces to linear interpoiatioQ 
alotnjg the edges of the unit square, and thus Is, mathematically reasonable. Further, 
this mdhod may be applied* tn thg thrMwiTnrit»ngimi'ft| ^ 

' hi Fig. ' 3 there is shown a unft cube interpolation' unit having ej^t vertices 
with oo-onUnxtesof 



10 



15 



20 



35 



iWi), (ii+l,yiX(), 

. .(xi^-+l,ft+l)/aod (xt+l^i+lA+1). • 

and including, a point P whia co-oxxliaatcs (xi+j%, yi+y^ lii-zf) at whkh the vahjc 
of U is to be interpolated. The cube is divided up into eight rectangular parakHo- 
pmedia by three planes which tnchide the . point P and are parallel to lu fe^^y j^e 
value U(3^^^) a the pdnt P is f<Rind by addii^ up tte 45 
each known value at each of the veztioes d the- unit, cube by the vi^ume of each* 
rectui^olar paraOdop^jcdon which b positioned opposite to that vertcx» thereby 
obtaioDag the foDowing formula (U): 
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+U?^-hl,yiA) X, (l-y,) (l^iA 
+U(xuy,+lA) (i-i,) 7t 
+U(»jriA+l) (1-xO (1-yr) 2, 

+U^ii+ljri+U) ir y, 
+U(*+l,yrhU+l) xr yi xc 

Agam^ this metfaod produces coaslsteot results at ^ vertfoes of the unit cube. 
Fbnlier^ along the edges of tbe unit cube k reduces to simple linear inteipoSarioo, 
and on the fac^ of tbe unk cube it reduces to the method of qoation (I). It'is furtiier 
5 clear that the value cAtained in the center of each face of tbe unit cube is the mean 

vahie of tbe known vahes at each vertex of that face, and the value obtained at the ^ 
center of the unit cube is the mean vahxe of the e4;ht known values at the vertices 
of the cube. Accordingly, dus merfiod ti seen to be mathematically zeastxiabte. 

However, this method has disadvantage! It lequires exg^ products to be 
10 fomae^ each of four vahies» apd addition therec^. Hence it is not always b^ for tn 
h^ speed calculation. 

There is another disadvantage in this method. Ahhough horn one mih cube 
to the next the huexpolated values ate comimums^ dietr derivative is not That ts» 
. the slope of the inteipolated values is discootinuoos from one unit cube to the next^ 

Le. the Ike of die inteipolated values bends sharply as we pass over the boundary* |c 
Thus in practice a sharp step of color values will be apparent in the finished picture, 
and the cubic structore of tibe .memory will show, to the detriment of qualinr. This 
effea can become quite serious. Fig. 4 shows a distnbution of the interpolated vahies 
obtained accordmg to the formula (I) which has a saddle forw, which shows the 
20 aforementioned inconvenience dearly. An even continnous line of interpolated vahxes 20 
in the unit stjuare AtBtCiDi is oboined, and also in die unit square AACfity but 
between these two square^ at tiieir common border, the derivative of the intcxpolated 
vahies is discontinuous. 

. The prior art incerpolation methodic and their disadvantages, have been explained 

^ above. A meduxl acootdhig to the present invention wiH now be described. 

In Ft^ St showing ti» twiHlknensional case^ two adjacent interpolation regions 
At^CtDs and AjBtQA are shown. The centezs of tfae^ unit squares ate designated 
by and 0«, and interpolated values at diese points are derived as averages of die 
functhn vahitt at die four comen of die square*. Tbm interpolation is conducted 
30 hneariy m «ch of die tnengles AAB„ BjO^Cu QCDi, DtO^Ai, A,OA, BAC 
QiOtDjt, and D,OtA«. That is, the point at whidi die value la to be interpolated is 
first checked to detnmme vrfuch of diese trian^ it hlh hxo, and then the value at 
ttie point is determmed by interpolation in die trkn^ in a fashion analogous to 
Fig. 2, by drawing lines &om ttsc point to the comers of die trian^c, and dicn 
calwlatmg die value of die function at die pomt as a wdghted sum of die values 15 
at comen of die trian^ givmg each value at a comer a wci^itiag of die lado 
<rf die area of a second triangle whose comers are the point and die odicr two coiners 
of die triangle, and die area of die trian^ In this mcdiod die magnitude of die 
discontinurty in die derivative of die inteipolated values firom one intcrpoUdon leaon 
to die next is much reduced, 

Now, considering Ae diree-dimcnsional case, die basic interpolaticm mediod in 
a tetrahcdtal ydumc will be explained widi reroea to Rg. 5. Let ABCD be a tctra- 

^'^i J'S'r.*^ ^'^^ " * P^"* ^ ^ <rf ^ function U to be 
mtwpolated is known. The value at point P, mtemal to die tetrahedron, is calculated 
as f Qflows :dr»wlm« from each vertexA.B,C,andD dmnigh die point Ptomeel * 
Sic' U(P)V ^n^^^A^n in A', C, amliy. Ttei Se interpolated 
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Nov the ratio of tfac volumes of the tetrahedra PBCD and ABCD, for example, 
is the same as the ratio of the heights of P and of A from the plane of BCD. 

Refentng to Fig. 7) theie is shown a unit cube interpolation vohune ABCDEFGH, 
and the vatoes of the function U of three variables are assumed to be known at the 
vcrtioes of . the cube. AB dnee varbtions of the piesou method dep^ upon dissecting 
this • cube into tetiaheldra w|iose veitioes aie either vertices of- the cube, centers of 
face* of the cub^ or the center of dte cube. Then a setiea of conqparisons are made 
to determine which of these tettahedra contains the point at which the value of the 
funcdon is re<jmied to be Interpolated. Once this is detemined, the value is then 
kterpolated within that tetrahedron according to die method described above^ using |a 
analyticai geometry; It will be lealiKd that it is ;nadiematicany teasooable to 
interpolate, initially, the values of the function at centers of faces of the cube as die 
aver age o f the vahies at the four comers of die faces, and the value of the function at 
die center of die cube as die average of the valtjes at aU eight vertices of the cube. 
Thus for each vertex of each tetrahedron of die dissection of Fig. 7 die value of die 15 
function b known, ^nd dierefiac die mcdtod aiustrated in Fig. 6 can be applied for 
mterpolatioiL' 

m die dteection of Fig. 7 diere exist twenty-four tetrahedra, of each of wUdk ■ 
ow vertex is die center of die cubic unit; two vertices are vertices of the cubic unit 
whoh are c onn e ct ed by an edge of die cubic uni^ and die fourth vertex is the cemer 
of 2 sqpMTc face of die cubk unit (me edge of which face is dte said edge of die cubic 
Mft The twenty-four tetrahedra are all isomorphic. One of diem is shown in Rg. g. . 

pe centen of die £aces of die cubic unit have be^i labelled as Q, O, and 

" ^- tetrahedron OABQ, Is iUustrated The plancs.OAB, . 
QiAB, OQxB, and have die equations y,-2f=(>, z,=0,.x,+y,-l=0, and 25 

j-y^^ xespecGvcIy. Thcicfore tt is dear that die condition for die pomt P to He - 
hmde dtf tecrahtttom OABQ,' is diat y,-2f>0, x,+y,<0. and x,^y,>0. (Of 
cotoe If >(i lqr deM^on). Provided diese con^^ 
vahie of die fancttra may be calaikifid as oudin^ 

U(A); X ratio of volumes of POQiB and OABQ, 
• ■ ' • • • ... 

UCB) X ratio of volumes of POQxA and OABQ. 

: . . • ^. . . . ...^ ...... 

. /^(Q^) ^ <^ o^.fPAB pAB 

UfO) X ratio of volumes of PQiAB aial OABQi 
. : . _ . =U(A) [i-x^-y,] +U(B) W-yJ 

^ . +m) Wyr-Xr] tU(0) (2i] (OT) .30 

This is became die ratio of the vohmies - of the abovcmentioned tetrai^^ 
, pointed out above, is the ratio of their heigfatst, and the etpiations of their faces axe 

as statcd.above. • ^ 

Sinoiilar results bold wheii the pohit P a in the odttr kterpoiation tetrahedra. 
^ A ioiniplete tabte of the' conditions for discrimination of which tetrahedron contains 35 
die point P, and of die fkcton iriuch are used for calculation of the interpolated value 
; in each case, is sEown as Table L Using diis table, by testing the oonditicns that 
are mat parendiesized, it is possible to characterise the tettah»lnm v^iiicfa contains 
the point P, and accordingly it is not necessary to test die pamdiesized conditions. 
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It IS readily ^mdemood diat the calculadon is far simpler in practice than tiie 
method of the abovoDentioned fomula (II). Further] in this method, the dis* 
comimitt^ across the bordecs between one unit cabe and the not are much leduoed, 
since the values near the hot of the tm^ cube axe much more daminatcd by the . 
5 values at the four comers of the hot than hi the pxior art method of (II). 5 
In fact the color pictuxe output tlffiaJs in the memory commonly vaxy motm- 
tmikallyj and tfietefore a more simple and coarse method ol fntapdatioD tfian the 
one oudhied above may well be sati^actozy m a particalar case* Therefore the method 
of Figs. 9 and 10 may wdi be acceptable^ ahhcmgh it b not quite so accnrate as the 
10 mediod of formula (III). In Fig. 9 is shown a dissection of the unit oibe mto iQ 
tetiahedta aQ of those vertices axe vertices of the unit cube. Thus this method ha;' 
the advantage that no avcragmg of vabes at the vertices of the unit cube a necessary 
h order to dfitffnnmit vahies at the centezs of the faces of the unit cube and at its 
center. 

15 Hie unit cube is dissected into sii tetraheda by three planes which have . a Ime 

hi commoo whidi Is the kmg diagonal of the imit cube^ and eadi plane is farttnwf 
to the other two at 50* and contams two edges and four vertices of the unit cube. A 
npical one of the sec tetrahedra is iOusixated ta Fig. la In ^ case die con^itiwig 
the point P to He within this t^rahedroa are that it>7t>Zu as can be easily 

20 worked out using solid geometry^ as before; In the same way the imerpolated value jf\ 
U(P)i8eqoalto 

U(A) [1-x,) + U(B) [i^yj + U(C) (yr-Zr] + U(D) £^ 

Similar diyw'm tnating conditions and calculating &ctors can be worked out for 
die other £ve tetrahedra. Table 2 shows ^ ooo^ete set It is readily ai^>reciited 
:25 that this variation of the method is easier h calculation than the method of fbnmila 25 
(HI), albeit at a slight loss in accuracy. 
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. of the method} iUustnted ^ Fig. 7. However this vaxiatiiHi will only be sliglit in the . 
case of monotcnic fimctfons, and it quite tolcrBt?le. 

WHAT WB CLAIM IS:— 

^- A Imear toteipdatmg mediod for cdor signals in a memory of a picture 

5 leprodudng maehioei compiismg stmiog ^ipiopnate values of ooior pictnic output ' c 
ngoals contsponding to ocrtaiD stepped values of color input signals in the memdiy' 
addressed in a tfaree-dxmeosional £uhk% and inteipolatiiig values of color output 
signab at points ^rfnch axe betsvten nid values by: 

dividmg up the cohk interpdatfon unit of the memoiy which is otmstttuted by 
10 a single step of each of the color input s^nals-ioto a phnnli^ of tetrahedra whose iq 
vertices ait either vertioes of the cubic uni^ oenten of its face% or its center point; v 

calculating the color output signal at each vertex of these tetrahedra whidi is a 
center of a hcc of the cubic uni^ if iny* by averagiatg the values of the color output 
signal at the four vertices which are comen of said face^ and at the center point of 
dK cubic unit by averaging ^ vahies ci the cohv output signal at all eidtt of the 15 
vertices of the cubic unit; 

dftfeTmrn fog vdiidi of these tetrahedra xndudes the interpolation point at i^iich 
the vahie of the color output signal is to be hnetpolatied; and 

deriving the hxtetpdated vahic at the interpolation pomt as a wef^to! sum of 
20 die vahia at die four vertices of the determined tetxahediQtt» the vahie at eadi vertex ^ 
being ghen a weight corresponding to the ratio of the vobme of t secosnd tetrahedron 
, whose vertices are the interpolation point and die other three vertices of the deter^ 
ffiined tetxahedrtm to that of the determined tEtrahednm. 

X A mediod as in Claim 1, whcreia the cubic unit is divided into twenty^fbor 
^ ' tetrahedra, of each of which one vertex is tiie center of the cubic unk, two vertices 
are vertices of the cubic unit which are connected by an edge of the cubk uiih; and 
die fourth vertex is the center of a square &oe of die cubic unit one edge ol which 
face is the said edge of the cubic unit 

3. A medux! as in Cairn 1, wherein die cubic unit is divided into six tetrahedra 
by three planes which have a line in common, said line being a Itmg Hw ghn a] of the 
aibic uni^ and each plane containing two edges and four vertioes df die cubic unit 

4. A mediod as In Claim 1, wherein die cubic unit is divided into five tetrahedra 
by four pknes, each of which contaks eiactiy three vertices of the cubic imil^ said 
planes mtexsecdng one anotiier aioQg Imes whkh are diagonals of faces of die cubic 
tmit 

^5. A linear interpolating mediod for o^our signals in a memory of a, picture 

rtproducing machine as datmed ia CUka I and substantially as hereinbefore described 
wiA reference to and as illustrated m FJgs. 5, ^ 7 tod 8, 9 and 10 or 11 and U of 
die accompanymg drawings. 

H: L. COTTRHLL & Ca„ 
Chartered Patent Agems, 
22 Whitefiiargate^ 
HaH HUl 2BX. 
Agents fbr the Applicants. 
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